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Abstract. We reanalyzed high quality spectra of 158 stars of spectral types AO-Fl and v sin; between 60 and ISOkms"'. 
Using a Least Squares Deconvolution technique we extracted high S /N broadening profiles and determined the loci of the 
Fourier transform zeros qi and q2 where the S /N-ralio was high enough. The values of v sin / were redetermined and found to 
be consistent with the values derived by Royer et al. (2002). For 78 stars q2 could be determined and the ratio ^2/^1 was used as 
a shape parameter sensitive for solar-like differential rotation (Equator faster than Pole). 74 of the 78 stars have values of ^2/91 
consistent with solid body rotation; in four of the 78 cases, values of q2/qi are not consistent with rigid rotation. Although 
these stars may be binaries, none of the profiles shows signatures of a companion. The Fourier transforms do not indicate any 
distortions and the broadening profiles can be considered due to single objects. One of those candidates may be an extremely 
rapid rotator seen pole-on, but for the other three stars of spectral types as early as A6, differential rotation seems to be the most 
plausible explanation for the peculiar profiles. 
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1. Introduction 

The substantial difference between photospheres of solar-type 
stars and A-type stars is the existence of a convective envelope. 
Due to the ionization of hydrogen the cooler late-type stars har- 
bour convective envelopes where turbulent motions of the pho- 
tospheric plasma can occur Stars of spectral types earlier than 
about F2 have no or only very thin convective envelopes and 
properties of granular flows change fundamentally. 

The generally accepted activity paradigm places the stel- 
lar dynamo believed to cause stellar activity at the bound- 
ary between the convective envelope and the radiative core. 
Differential rotation drives the dynamo action by winding up 
and amplifying the magnetic flux tubes. The interaction of 
magnetic fields, differential rotation and the convective en- 
velopes are believed to be ultimately responsible for stellar ac- 
tivity. 

In a series of publications an onset of convection was 
searched for. It is generally accepted that the onset of stellar 
activity occurs between spectral types A7 and F5 depending on 
the observational strategy. Wolff et al. (1986) studied Cii and 
He I emission and placed the onset of activity near B-V - 0.28, 
i.e. around spectral type FO. Schmitt (1997) concluded from 



X-ray data that coronal emission is universal in the spectral 
range A7 to G9 implying an onset of activity around spec- 
tral type A7. Hotter stars are expected to harbour shallow con- 
vective envelopes, these stars have higher convective velocities 
which peak at about A3 until convection disappears altogether 
at about Al (Renzini et al. 1977). Gray & Nagel (1989) directly 
searched for the onset of convection analyzing line bisectors of 
slowly rotating stars. In their targets the Doppler-shift distri- 
bution of the granulation dominates the broadening of spectral 
lines and a bisector reversal was found around spectral type FO. 
Stronger asymmetries were found in the stars at the hot side of 
the boundary indicating higher convective velocities. 

Although stellar activity is not observed in early A-type 
stars, it is not clear whether differential rotation may take place 
in early-type stars. The absence of activity may simply reflect 
inefficient coupling of surface magnetic fields and the lack- 
ing interface between the radiative core and a convective enve- 
lope. There is no reason to believe that rapidly rotating A-stars 
should rotate rigidly. In case of the late-type Sun we know that 
the surface rotation law can be approximated by 



= QEquator(l - ff sin^ 0, 



(1) 
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* Based on observations collected at the European Southern 
Observatory, La Silla 



with / the latitude and a© ~ 0.2 as derived from Sun spots. 
In the solar case the Equator rotates about 20% faster than the 
Poles. Gray (1977) searched for differential rotation in line pro- 
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files of six A-stars finding no indications of differential rotation 
within his error bars. 

Also using line profiles, Reiners & Schmitt (2003) found 
signatures of differential rotation in a sample of F-type stars. 
The earUest object in their sample indicating differential rota- 
tion is of spectral type FOIV/V. Applying the method used by 
Reiners & Schmitt (2003), we search for signatures of differen- 
tial rotation in a large sample of A-star spectra. The results are 
presented in the following. 

2. Observations and data analysis 

The spectra were observed with the ECHELEC spectrograph 
(ESO/La Silla) and are part of a larger sample collected in 
the framework of an ESO Key Programme. These observa- 
tions were aimed at the determination of fundamental parame- 
ters of early-type stars observed by HTPPARCOS (Gerbaldi & 
Mayor 1989). The total sample is described by Grenier et al. 
(1999) who measured radial velocities, and Royer et al. (2002, 
hereafter RGFG) who derived rotational velocities from these 
spectra. The observed spectral range spans from 4210 to 4500 
A. The hnear dispersion is about 3.1 Amm the slit width 
of 320 fxm corresponds to V.'52 on the sky, and the resolving 
power is about 28 000. 

To search for the spectral signatures of stellar rotation 
laws, broadening profiles were derived by applying a Least 
Squares Deconvolution process (LSD). After constructing a 
(5-template comprising the strongest 150 lines taken from the 
Vienna Atomic Line Database (Kupka et al. 1999) and ac- 
cording to stellar temperature, a first-guess broadening profile 
was deconvolved using each pixel as a free parameter in the 
fit. Since theoretical line depths match the observational ones 
poorly, the equivalent widths of the incorporated lines were 
optimized in a second step while leaving the broadening pro- 
file fixed. During a few iterations the broadening profile and 
the equivalent widths were optimized. Using this technique the 
spectral lines are effectively deblended, the information con- 
tained in every spectral line is used and the signal-to-noise ratio 
is significantly enhanced. Consistency of the fit is checked by 
comparing theoretical line depths to the derived ones. 

Following Reiners & Schmitt (2002) we Fourier trans- 
formed the broadening functions and measured the position of 
the first and second zeros {q\,q2)- The ratio q2/q\ is a robust 
observable for the shape of a rotational broadening function 
and a direct indicator for solar-like differential rotation with 
the Equator rotating faster than the Pole (cp. Reiners & Schmitt 
2002). We measured the ratio q2/q\ for all stars the LSD pro- 
cedure yielded a stable and symmetric broadening function for. 

The spectral quahty used in this analysis in principle was 
sufficient to follow the Fourier transformed broadening func- 
tions to the second zero q2 in stars with projected rotational 
velocities in the range 60 km s~' < v sinj < 150kms~^. 

3. Results 

The outhned method was apphed to the spectra of 158 stars 
of spectral types AO-Fl. The rotational velocity v sin i was de- 
rived from the first zero qi. For 78 of our sample stars the ratio 
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Fig. 1. Comparison of the v sin ; values with the ones derived by 
RGFG. The solid line is the one-to-one relation, and the linear re- 
gression between both data sets is given in Eq. 2. 

q2/q\ could be determined. For the disregarded 80 stars, either 
data quality was insufficient or the derived broadening function 
showed obvious peculiarities probably due to binarity. 

The V sin / values are compared with the results from RGFG 
in Fig. 1, for the subsample of 78 stars. RGFG also determined 
their values of v sin; from qi but used selected absorption 
lines while an "overall" broadening profile was deconvolved 
here. The consistency of both scales is robustly estimated using 
GaussFit (Jefferys et al. 1998a,b). The resulting hnear relation 
is 

V sin ; - 0.99±o.05 v sin /rgfg - 1.6±4.6, (2) 

both scales are in good agreement. The systematic differences 
occuring due to the different ways obtaining the broadening 
profiles can be neglected in our sample and should not be dis- 
cussed here. 

We measured the second zeros of the Fourier transformed 
broadening profiles to calculate the ratios q2lq\ - The results are 
plottedinFig. 2, typical errors are of the order of A^2/^i ~ 0.1. 
A rigid rotator is expected to yield a value of q2lqi between 
1 .72 and 1 .83 assuming a linear limb darkening law (indicated 
by dashed lines in Fig. 2). For the stars of our sample lin- 
ear hmb darkening coefficients between 0.5 and 0.75 are ex- 
pected during their time on the Main Sequence (Claret 1998). 
Assuming a limb darkening parameter of e = 0.6 rigid rota- 
tion would yield q2lq\ = 1.76 (solid line in Fig. 2). The results 
that are consistent with a value oi q2/q\ = 1.76 within the er- 
ror bars are indicated by small crosses in Fig. 2, for the sake 
of readabiUty no errors are plotted for them. The second zero 
q2 can only be determined in spectra where the signal exceeds 
the noise level beyond q2, i.e. when a second sidelobe is de- 
tectable. In our case the ampUtude of the second sidelobe is at 
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Fig. 2. Derived values of qi/qi plotted against v sin (' as derived from 

the first zero of the Fourier transform (with l-cr uncertainties). The re- 
gion between dashed lines is consistent with solid body rotation for 
arbitrary Umb darkening. For linear limb darkening with e = 0.6, 
gi/qi = 1-76 is expected (solid line). Crosses indicate results con- 
sistent with solid body rotation, typical errors are of the order of 
^<iiUl\ =0.1 (not plotted). Four stars not consistent with (/2/i7i = 1.76 
are indicated by solid circles, error bars are plotted for them. 



the noise level for many stars, and these measurements of q2 
must be interpreted as lower limits; thus some of the measure- 
ments of q2lq\ plotted as crosses in Fig. 2 are essentially lower 
hmits. For 74 of the 78 stars analyzed the broadening profiles 
are consistent with solid body rotation. 

Four of our measurements are not consistent with qilqi - 
1.76, they are indicated by full circles in Fig. 2 and errors bars 
are plotted for them. For three of those — HD 6869 (A9V), 
HD 60555 (A6V) and HD 109238 (FOIVA^) — the values 
of qilqi are significantly smaller than 1.7. The fourth star 
(HD 44892, A9/F0IV) has a value of qilqx marginally con- 
sistent with qilqx > 1.7 within its error bars. We will discuss 
this star in Sect. 3.2. 

The Fourier transforms of HD 6869, HD 60555 and 
HD 109238 are plotted with error bars in Fig. 3. Overplotted 
are the Fourier transformed line profiles of stars with sim- 
ilar values of v sin? that are consistent with rigid rotation 
(li/qi - 1.76). While different velocity fields, e.g., turbulence, 
may influence the amphtudes of the sidelobes, the zeros of the 
Fourier transform arise from rotational broadening (cp. Gray 
1976). One mechanism known to change the ratio q2/q\ in the 
manner found in HD 6869, HD 60 555 and HD 109 238 is solar- 
Uke differential rotation. The strength of differential rotation in 
terms of the parameter a in Eq. 1 can be calculated from qi/qi 
(Reiners & Schmitt 2003), and the respective values of a are 
given in Table 1 together with the spectral types and v sin i of 
the four suspected differential rotators. 
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Fig. 3. Fourier transforms of HD 6869 (top panel), HD 60555 (mid- 
dle), and HD 109238 (bottom panel) plotted with error bars. These 
stars show extremely small values of q2/q[. The spectra of HD 29 920 
and HD 145 876 have values of qi/qi as expected for rigid rotators and 
are plotted with dashed lines for comparison in the top (HD 29 902), 
middle and bottom (both HD 145 876) panels, respectively. 



3.1. Binarity 



The spectra we used were also studied in the aim of deriving 
radial velocities by Grenier et al. (1999). The authors gave 



a "shape" flag for the correlation functions for our four 
candidates: 
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Table 1. Derived values of v sin; and (^2/92 for the four stars with values of ^2/91 significantly smaller than 1.76. Also given are strength of 
differential rotation in terms of a (cp. Eq. 1) and required values of equatorial velocities v'e_,igid and inclination angles ; if the value of ^2/^1 
is explained by rapid rotation seen pole-on (Sect. 3.2). Distances in pc, ROSAT X-ray luminosities and v sin i derived in RGFG are given in 
columns 8-10 



HD 


Type 


V sin! 


?2/?l 


a 


Ve,rigid 


i 


d 




V sini 






[kms->] 






[kms-i] 




[pc] 


[W m-2] 


(RGFG) 


6869 


A9V 


100 ±6 


1.55 ±0.08 


0.28 + 0.10 


(460) 


(13°) 


87 


570 


121 


60 555 


A6V 


115 ±7 


1.54 + 0.07 


0.29 + 0.08 


(470) 


(14°) 


134 




114 


109 238 


FOIVA' 


103 ±4 


1.51+0.11 


0.32 + 0.13 


(500) 


(13°) 


133 




125 


44892 


A9/F0IV 


152 + 5 


1.64 + 0.10 


0.16 + 0.16 


400 


22° 


160 




145 



- HD 6869: "probable double", 

- HD 44 892: "probable double", 

- HD 60 555 : "suspected double", 

- HD 109 238: "probable double". 

Using our deconvolution method, however, we find indication 
for double peaks in the broadening functions of 23 of our tar- 
gets, but not in the spectra of these four stars. Note that this 
does not mean that the 135 others are single stars since the 
luminous A-type stars dominate spectra of, e.g., binaries con- 
sisting of an A-type and a G-type star. The G-type spectrum 
will easily be hidden in the light of the A-type star. To be com- 
plete we checked the shape of the correlation function by cross- 
correlating our template with the spectra and found no indica- 
tions for binarity either. 

For HD 44892 and HD 109 238, the hterature gives hints 
about a single star status. There is no evidence of binarity for 
HD 44 892, neither in HIPPARCOS data nor in Speckle obser- 
vations (Mason et al. 2001). It can be considered as single with 
a high level of confidence, and its spectrum is surely not af- 
fected by any significant contamination. HD 109 238 is part of 
the sample observed by Abt & Morrell (1995). The MK classi- 
fication they derive for this object is FOV, with no suspicion of 
spectroscopic binarity. 

For the two other objects, HD 6869 and HD 60555 two 
spectra are available in our data set. Individual observations of 
HD 6869 and HD 60 555 are separated by 383 d and 767 d, re- 
spectively. Inspection of the broadening functions derived from 
the individual spectra yields no indication of variability due to 
relative motions of binary components. Coadded spectra were 
used to derive the values of ^2/^1 for both stars. Both targets 
are indicated as binary stars in the literature and are discussed 
below. 

3.1.1. HD 6869 

HD 6869 is a binary system seen with an angular separation 
of about 1'.'2 (ESA 1997; Horch et al. 1996). With such a low 
separation, compared to the width of the sUt, the collected 
spectra are expected to be contaminated by the light of the 
secondary. The individual magnitudes in HIPPARCOS and 
TYCHO bands are given by Fabricius & Makarov (2000) for 
both components and are respectively: 



- for the primary: Hp = 7.190 + 0.003 mag, 

Bt = 7.42 ± 0.01 mag, 
Vt =7.16 ±0.01 mag, 

- for the secondary: Hp = 8.790 + 0.014 mag, 

Bt = 9.34 + 0.01 mag, 
Vt = 8.76 + 0.01 mag. 

No significant interstellar absorption is expected (Lucke 1978) 
from the position of the star in the solar neighborhood (I - 
294°, b = -70°, d = 87 pc). Considering no reddening, 
the temperatures of each component, derived from the Bt - 
Vt color index (Appendix A), are respectively 7500 K and 
6000 K, according to Eq. A.l. The magnitude difference in Bt 
corresponds to a flux ratio of about 5.5 around 4450 A, tak- 
ing into account the effective temperature of both components. 
Given the angular separation and the large magnitude differ- 
ence, the spectrum of HD 6869 can be considered as a single 
star spectrum of the dominating A-type star. 

3.1.2. HD 60555 

This star has been very little studied. Its spectral type in the cat- 
alogue published by Houk (1982) indicates a composite spec- 
trum: A5/7V-i-(F). Grenier et al. (1999) flagged this star as vari- 
able in radial velocities from these two spectra. The ratio E/I of 
external to internal error on radial velocity is 2.93. The absorp- 
tion spectrum of HD 60 555 in the observations does not show 
any evidence of multiplicity, especially in the Fourier transform 
no indication of contamination due to the light of a secondary 
is apparent. 

3.2. Extremely fast rotation 

Alternative to differential rotation the shape of the broadening 
function and the value of qi/qi can also be affected by very 
rapid rotation and gravity darkening possibly observed pole-on 
(cp. Reiners 2003). Flux is redistributed from the Une's wings 
to the center when the Equator becomes cooler due to gravity 
darkening. As far as the lines considered are not dominated by 
temperature and gravity variations over the stellar surface — 
which is the case, e.g, in the weak lines of early A-type stars 
as shown by Gulliver et al. (1994) — the value of q2/qi is di- 
minished by this effect. According to Reiners (2003), the ratio 
qi/qi then only depends on the equatorial velocity Ve and on 
the gravity darkening law. We assume a linear gravity darken- 
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ing law according to Claret (1998) and calculate the equato- 
rial velocities Ve.rigid required to produce the measured values 
of q2lq\ assuming solid body rotation for the four suspected 
differential rotators. The results and the respective inclination 
angles i are given in columns six and seven of Table 1 . 

For HD 6869, HD 60555 and HD 109238 the values of 
Ve,rigid are larger than breakup velocity; for these stars rapid 
soUd body rotation can be ruled out as the mechanism solely re- 
sponsible for the diminished ratio q2/q\- In case of HD 44 892 
the rotational velocity required for the measured ratio q2/qi = 
1.64 is of the order of breakup velocity. Thus differential rota- 
tion as well as rapid solid body rotation are the two possible 
explanations for the measured profile shape of HD 44 892. 

4. Conclusions 

We reanalyzed high quahty data formerly discussed by Grenier 
et al. (1999) and Royer et al. (2002), with the aim of search- 
ing for differential rotation in early-type stars. With an iterative 
Least Squares Deconvolution method we obtained high qual- 
ity broadening profiles of 158 stars with projected rotational 
velocities in the range 60kms"' < v sin/ < 150kms"'. We 
disregarded the profiles of 80 of them due to obvious asymme- 
tries or multipUcity. For 78 stars the broadening profiles appar- 
ently reflect the rotational broadening law. Profile distortions 
were analyzed in terms of the ratio of the first two zeros of 
the Fourier transform qilqx. Within the errors, 74 of the 78 
measured profiles are consistent with the assumption of rigid 
rotation. Due to data quahty many measurements must be con- 
sidered lower limits and from this sample no conclusion can be 
drawn concerning values of q2lq\ possibly larger than 1.8. 

Unfortunately, for none of our sample stars interferometric 
measurements are available to our knowledge, e.g., Altair and 
Sirius have rotational velocities outside our range determined 
by data quahty, and they are not contained in our sample. 

Four stars are analyzed in detail, the profile of the A9/F0IV 
star HD 44 892 is only marginally consistent with rigid rota- 
tion. It is hkely that its profile is distorted either by differential 
rotation or by very rapid rotation seen pole-on; in the latter case 
HD 44 892 would be the first star that directly shows signatures 
of gravity darkening in mean profile broadening as proposed by 
Reiners (2003). Comparison to interferometric results would be 
interesting especially for this star. 

The broadening functions of the three stars HD 6869 
(A9V), HD 60555 (A6V) and HD 109238 (FOIV/V) are not 
consistent with rigid — even very rapid — rotation since their 
equatorial velocities would be larger than breakup velocity. 
Although some authors suspect these stars being binaries, in 
our high quality spectra we find no indications of multiplicity 
neither in data nor in Fourier space. Since contamination due 
to secondaries are easily visible in Fourier space — where no 
sharp zeros should occur in case of the profile being a sum of 
two — we consider the spectra single star spectra. Differential 
rotation seems to be the most plausible explanation for the ob- 
served profile distortions. For these three stars the Equator is 
rotating about 30% faster than the polar regions. Thus we con- 
clude that significant differential rotation seems to take place 



even in early-type stars not harbouring deep convection zones, 
the earliest object is the A6 dwarf HD 60 555. 

If differential rotation is the driving mechanism for stel- 
lar activity, these stars should be active, too. X-ray emission 
from HD 6869 was detected with the ROSAT mission and the 
other stars may also be X-ray sources but were simply too far 
away for a detection. Whether differential rotation is a common 
phenomenon in these stars cannot be answered by this work 
since only very strong differential rotation is detectable with 
our method. The finding of strong differential rotation among 
A-type stars indicates that there is no abrupt change in rota- 
tional laws of stars around the boundary where surface convec- 
tion sets in. 

Acknowledgements. A.R. acknowledges financial support from 
Deutsche Forschungsgemeinschaft DFG-SCHM 1032/10-1. This 
work was partly supported by the Swiss National Science Foundation. 

Appendix A: Bj - Vj versus Teir 

The effective temperature can be approximated using the color 
index Bj-Vt- This calibration is built from the catalogue given 
by Cayrel de Strobel et al. (1997). They compile [Fe/H] and 
Teff determinations from the literature. These data together with 
TYCHO colors for these stars allow the calibration of Teff as a 
function of Bj - Vj. The resulting formula is: 

log(reff) = 4.0132- 0.87573 X -I- 2.2194 (A.l) 

-5.0087 + e.ieie / - 5.1069 

-1-2.0638 / - 0.39629 x' + 0.024548 x^ 
with x = Bt - Vt- 
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